Drag and diffusion coefficients of B mesons in hot hadronic matter 
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The drag and diffusion coefficients of a hot hadronic medium consisting of pions, kaons and eta 
using open beauty mesons as a probe have been evaluated. The interaction of the probe with the 
hadronic matter has been treated in the framework of chiral perturbation theory. It is observed 
that the magnitude of both the transport coefficients are significant, indicating substantial amount 
of interaction of the heavy mesons with the thermal bath. The results may have significant impact 
on the experimental observables like the suppression of single electron spectra originating from the 
decays of heavy mesons produced in nuclear collisions at RHIC and LHC energies. 
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I. INTRODUCTION 

The suppression of the transverse momentum 
(pt) distribution of hadrons produced in nucleus- 
nucleus relative to (binary scaled) proton-proton in- 
teractions at the Relativistic Heavy Ion Collider 
(RHIC) [l| has been used as a tool to understand the 
properties of matter formed in such collisions. The 
large value of the elliptic flow of hadrons measured 
at RHIC along with the suppression of the high px 
hadrons mentioned above indicate that the matter 
might have been formed in the partonic phase with 
liquid like properties characterized by low value of 
shear viscosity (ry) to entropy density (s) ratio, ry/s 
with a lower bound of rj/s ^ 1/Att 

In addition to elliptic flow (U2) and nuclear sup- 
pression (Raa) of light hadrons these quantities have 
also been measured for the single electron spectra 
originating from the decays of the open charm and 
beauty mesons produced at RHIC coUisions [1, 
The advantages with heavy mesons are two-fold. 
Firstly, they contain either a charm or a beauty 
quark which is produced very early and hence can 
witness the evolution of the partonic matter since its 
inception until it reverts to hadronic matter through 
phase transition/cross over Q and secondly, the 
heavy quarks do not decide the bulk properties of 
the latter. Therefore, charm and beauty quarks are 
considered to be efficient probes for the characteri- 
zation of the partonic phase. In most of the earlier 
works (cl-fTsI aimed at extracting the properties of 
quark gluon plasma (QGP) by analyzing the Raa 
and V2 of heavy flavours the role of the hadronic 
matter was ignored. However, for the character- 
ization of QGP the interactions of heavy flavours 
with hadronic matter should be taken into consider- 
ation and the effects of hadrons must be subtracted 
out from the observables. Though a large amount 
of work has been done on the diffusion of heavy 
quarks in the QGP the diffusion of heavy mesons 
in hadronic matter has received much less attention 
so far. Recently the diffusion coefficient of D meson 



has been calculated using heavy meson chiral per- 
turbation theory jioj and also by using the empirical 
elastic scattering amplitudes [20] of D mesons with 
thermal hadrons. The Z)-hadron interactions also 
have been evaluated using Born amplitudes 12 ll and 
unitarized chiral effective Dn interactions |22l ). It 
has been found that the contributions of B meson to 
the single electron spectra dominate over those from 
D meson for large transverse momentum, px > 5 
GeV [2^ (see also [13]). Moreover, the future experi- 
ments are progressing toward precision measurement 
over a wide range of kinematical variables. In view 
of this the use of B meson as a probe to extract the 
properties of matter at high temperature assumes 
importance. 

In the next section we discuss the formalism 
adopted to evaluate the drag and diffusion coeffi- 
cients of the heavy flavoured mesons in a hadronic 
matter consists of pions, kaons and eta. Results are 
presented in section HI and section IV is dedicated 
to summary and discussions. 



II. FORMALISM 

In the present work the drag and diffusion co- 
efficients of the B meson propagating through a 
hot hadronic matter are evaluated within the am- 
bit of Heavy Meson Chiral Perturbation Theory 
(HMxPT) in LO, NLO and NNLO approximations. 
We also revisit the transport coefficients of D me- 
son in a similar theoretical framework. We consider 
the elastic interaction of the B meson with thermal 
pions, kaons and eta in the temperature (T) range 
100 — 170 MeV. Detailed analysis of the experimen- 
tal data on the hadronic yield in heavy ion collisions 
show that the value of the temperature for the chem- 
ical freeze-out of the system produced at RHIC en- 
ergies is about 170 MeV (see [25[ for a review). This 
indicates that the inelastic interactions which are re- 
sponsible for the change in the number of hadrons 
become rarer for T below ~ 170 MeV. Thus the con- 
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tributions of the inelastic collisions in evaluating the 
drag and diffusion coefficients of the hadronic mat- 
ter probed by the heavy flavoured mesons can be 
ignored. 

The drag (7) and diffusion (Bq) coefficients of the 
heavy mesons are evaluated using elastic interaction 
with the thermal hadrons. For the (generic) process, 
B{p) + h{q) — B{p') + h{q') [h stands for pion, kaon 
and eta), the drag 7 can be calculated by using the 
following expression [26| : 



7 = ViMlv^ 



(1) 



where Ai is given by 
f ( d^q 
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where the heavy-light pseudoscalar meson triplet 
P = (_B", i3+), heavy-light vector meson triplet 
P*, = {Bl^ ^ B*+ , B*+) and (...) denotes trace in fla- 
vor space. The covariant derivatives are defined as 
V^^Pa = df^Pa - PbTp and V^P^ = S^P] 
with a, b are the SU (3) flavor indices. 

The vector and axial- vector currents are respec 
tively given by = ^{u^dij^u + udf^u^) and = 



pt 

^ ab^b 



i{u^dfj,u — udf^u^) where u = exp(^). The unitary 



matrix $ collects the Goldstone boson fields and is 
given by 

JL 

_iLl + JL 
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To lowest order in $ the vector and axial-vector 
currents are: 



(6) 



gB being the statistical degeneracy of the B meson. 
The factor /(g) denotes the thermal phase space 
factor for the particle in the incident channel and 
l-|-/((jf') is the Bose enhanced final state phase space 
factor. From Eq. [2] it is clear that the drag coef- 
ficient is a measure of the thermal average of the 
momentum transfer, p — p', weighted by the interac- 
tion through the square of the invariant amplitude, 

The diffusion coefficient, Bq can be defined as: 



{{{p-p'?)) 



(3) 



Both the drag and diffusion coefficients can be 
evaluated from a single expression: 
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with an appropriate choice of T{pf). In Eq. |4] 
X{x, y, z) = + y"^ + — 2xy — 2yz — 2zx, is the 
triangular function. 

We start our discussion on the determination of 
the scattering amplitudes with the Lagrangian of 
Covariant Chiral Perturbation Theory (CxPT) in- 
volving the heavy B (or D) mesons [27] given by 



CcxPT = {V^PV^P^)~ml{PP^) 

-{D^,P*''Vt'P*'') +ml4P*''P*'') 



+ig{PXP^ - Pu^'P;^) 



(5) 



From the first term (or kinetic part of the P- 
fields) of CcxPT I the matrix elements for contact 
diagram in terms of Mandelstam variables (s, i, u) 
are obtained as 



Mb+t,+ = -Mb+^- = -7t;5(s - "") 
Mb+^o = MB + r, = Mb+k+ = Mb+k- = 



■[s-u) (7) 



These can be represented in the isospin basis as 
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where the isospin of the 5$ system appears in the 
superscript. Denoting the threshold matrix elements 
by T, these are obtained from ^ and are given by 
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One can reproduce these T-matrix elements in the 
isospin basis using the lowest order HMxPT La- 
grangian for heavy mesons containing a heavy quark 
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and a light antiquark of flavor a as given below 



trDiH^ri'ufH^) + .. .{10) 



where ff^^ 



2 V-" an / 



iPa"f^) and HQ 



-iPj^7^)ii^ and trjj denotes trace in Dirac 



space. In this formalism, since the factor ^/mp and 
^mp' have been absorbed into the Pa and P^^ 

fields, the threshold T-matrix element {T^) now 
has the dimension of scattering length ap whereas 
in C'xPT, we get a dimensionless T-matrix element 
(Tj^*). The relation between these two T-matrix 
elements and the scattering length ap is given by 



mpT^^ — 87r(m$ + nip)ap (11) 



The square of the isospin averaged T-matrix ele- 
ment is given by 



^|Tb$P- |Ts.|2 + |Tbk|2 



(12) 



where |Tb^|2 = 
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III. RESULTS 

We evaluate the drag coefficients of the _B-meson 
by using the momentum dependent and momentum 
independent matrix elements given by Eqs. H] and 
|9] respectively. The results are shown by the dot 
dashed and dashed lines in Fig. I. Inspired by the 
fact that the results for the two scenarios are not 
drastically different in the LO we proceed to evalu- 
at e the d rag c oeffic ient of heavy mesons by replacing 
X; |M|2 by J2 \Tf in NLO and NNLO also where the 
T-matrix elements will be obtained from the scatter- 
ing lengths. 

Liu et al [1^ have obtained the _B$ scattering 
lengths (see also [sO]) up to NNLO in HMxPT by 
using the coupling constant from recent unquenched 
lattice resuhs Using these NLO and NNLO 

results we estimate the isospin averaged drag coef- 
ficients of B mesons. The results are depicted in 
Fig. [T] The drag coefficient evaluated with NNLO 
matrix elements increases by 22% compared to the 
NLO resuh at T = 170 MeV. 

We now focus on the temperature dependence of 
the drag coefficient of P-mesons as shown in Fig. [TJ 
As mentioned before, 7 is the thermal average of 
the square of the momentum exchanged between the 
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FIG. 1: The variation of drag coefficients with tempera- 
ture due to the interaction of B mesons (of momentum 
100 MeV) with thermal pions, kaons and eta. The dot- 
dashed (dashed) line indicates the results obtained by 
using the matrix elements of Eq. [8] (T matrix of Eq. [9|. 
Solid (dotted) line indicates the results for NNLO (NLO) 
contribution 12911. 



heavy mesons and the bath particle weighted by the 
interaction strength through the invariant amplitude 
of the process. Therefore, with the increase in tem- 
perature of the thermal bath the kinetic energy of 
the hadrons increases. Hence the hadrons gain the 
ability to transfer larger momentum during their in- 
teraction with the B mesons - resulting in the in- 
crease of the drag coefficient. This tendency is ob- 
served in Fig. [T] quite clearly. The increase of drag 
with temperature is characteristic of a gaseous sys- 
tem. In case of a liquid the drag diminishes with 
T (except for very few cases). In this case a signifi- 
cant part of the thermal energy goes into making the 
attraction between the interacting particles weaker 
and once this happens the constituents move more 
freely resulting in a smaller drag force. Therefore, 
the variation of the drag with temperature can be 
used to understand the nature of interaction of the 
fluid. 

Since the diffusion coefficient involves the square 
of the momentum transfer it is also expected to in- 
crease with T. This is seen in Fig. [5J The drag 
and the diffusion coefficients are related through the 
Einstein relation as: 



Bo = MbiT. 



(13) 



where Mb is the mass of the B-meson. The tempera- 
ture dependence of the diffusion coefficient evaluated 
by using Eqs. |3] and the [T2] are displayed in Fig. [5] 
The difference between the results obtained from 
Eq. 0] and the Einstein's relation is about 6 — 7% at 
T = 170 MeV for the B meson momentum, p = 100 
MeV. This small difference illustrates the validity of 
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FIG. 2: Variation of diffusion co-efffcient as a function 
of temperature. The solid line indicates the variation of 
the diffusion coefficient with temperature obtained from 
Eqs. [3] and |4l The momentum of the B meson is taken 
as 100 MeV. The dashed line stands for the diffusion 
coefficient obtained from the Einstein relation (Eg. I13p . 



0.03 



0.02 



0.01 



0.03 



/ - 0.02 




0.01 



0.1 0.12 0.14 0.16 0.1 0.12 0.14 0.16 

T (GeV) T (GeV) 

FIG. 3: The variation of drag coefficients of D mesons 
with temperature due to interaction with thermal pions, 
kaons and eta in LO, NLO and NNLO approximations 
for interactions of D with thermal hadrons taken from 
Ref. 113] (left panel) and ^ (right panel). 



the Einstein relation in the low momentum (non- 
relativistic) domain. 

The energy loss oi a B meson moving through a 
hadronic system may be estimated from the relation 



dE 
dx 



IP 



(14) 



The magnitude of 7 obtained in the present calcula- 
tion reveals that the B mesons dissipate significant 
amount of energy in the medium. This might have 



crucial consequences on quantities such as the nu- 
clear suppression factor of single electrons originat- 
ing from the decays of heavy mesons. 

We also evaluate the D meson drag and diffusion 
coefficients using the interactions of D mesons with 
thermal hadrons discussed in Refs. \27] and 29] in 
LO, NLO and NNLO approximations. The results 
are displayed in Fig. [S] In the LO approximation 
the drag is similar for both the cases. However, for 
NLO and NNLO, the drag co-efficient evaluated us- 
ing the T-matrix elements obtained from the scat- 
tering lengths of Ref. |29] is slightly higher than that 
obtained from Ref. . 

The drag of D mesons in hot hadronic matter has 
recently been studied by using different approaches. 
While empirical scattering cross sections were used 
in Ref. [20j , the authors of Ref. [l^ used unitarized 
chiral effective Dtt interactions to evaluate the drag. 
We observe that the magnitude of the drag of D 
meson obtained in the present work is similar to that 
obtained in Refs. [l^l and [20|. The smaller value in 
the present case is due to the lower values of the D 
meson-hadron cross sections. 



IV. SUMMARY AND DISCUSSIONS 

In summary we have evaluated the drag and dif- 
fusion coefficients of open beauty mesons interact- 
ing with a hadronic background composed of pions, 
kaons and eta. It is found that the values of both 
the transport coefficients increase with temperature. 
The magnitude of the drag coefficient of the B me- 
son indicates that while evaluating the suppression 
of the high px single electrons originating from the 
decays of B mesons the effects of hadrons should be 
taken into account. Within the same formalism, the 
transport coefficients of the D meson has been cal- 
culated. The D meson drag coefficient is found to 
be lower than the values obtained in Refs. [20l. [2^. 

Some comments on the effects of the exclusions 
of inelastic channels and non-perturbative processes 
on the drag and diffusion coefficients are in order 
here. The inelastic channels do contribute to the 
scattering matrix through coupled channels via loops 
in the unitarisation procedure. However, in a ther- 
mal background such as in a heavy ion collision the 
scattering amplitude is weighted by phase space fac- 
tors which essentially control the rate of reactions. 
The fact that heavy ion collisions undergo chemical 
freeze-out at about 170 MeV means that number 
changing reactions will certainly be inhibited be- 
low this temperature. The presence of resonances 
makes the evaluation of scattering amplitudes a non- 
perturbative problem. Unitarisation of scattering 
amplitudes, preferably with loops evaluated with 
thermal propagators will certainly improve the reli- 
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ability of our results at higher energies. However, in 
the absence of any information about B* mesons so 
far, this exercise will however be far less constrained 
than the charm sector where the masses and widths 
of the excited states are known. 



Acknowledgment : 



SKD and JA are partially supported by DAE- 
BRNS project Sanction No. 2005/21/5-BRNS/2455. 



[1] I. Arsene et al. (BRAHMS Collaboration), Nucl. 

Phys. A 757, 1 (2005); B. B. Back et al. (PHO- [16 
BOS Collaboration), Nucl. Phys. A 757, 28 (2005); 
J. Adams et al. (STAR Collaboration), Nucl. Phys. [17 
A 757, 102 (2005); K. Adcox et al. (PHENIX Col- 
laboration), Nucl. Phys. A 757, 184,(2005). [18 

[2] P. Kovtun, D. T. Son and O. A. Starinets, Phys. 

Rev. Lett. 94, 111601 (2005). [19 

[3] S. S. Adler et al. (PHENIX Collaboration), Phys. 

Rev. Lett. 96, 032301 (2006); B. L Abeleb et al. [20 
(STAR Collaboration), Phys. Rev. Lett. 98, 192301 
(2007). [21 

[4] S. S. Adler et al. (PHENIX Collaboration), Phys. 

Rev. Lett. 98, 172301 (2007). [22 

[5] S. Borsanyi et al., JHEP 1011 (2010) 077. 

[6] H. van Hees, R. Rapp, Phys. Rev. C,71, 034907 [23 
(2005). 

[7] H. van Hees, M. Mannarelli, V. Greco and R. Rapp, [24 

Phys. Rev. Lett. 100, 192301 (2008). 
[8] C. M. Ko and W. Liu, Nucl. Phys. A 783, 23c [25 

(2007). 

[9] A. Adil and I. Vitev, Phys. Lett. B 649, 139 (2007). 
[10] P. B. Gossiaux and J. Aichelin, Phys. Rev. C 78, [26 

014904 (2008). [27 
[11] S. Wicks, W. Horowitz, M. Djordjevic and M. Gyu- 

lassy Nucl. Phys. A 784, 426 (2007). ' [28 

[12] S. K Das, J. Alam and P. Mohanty, Phys. Rev. C 

80, 054916 (2009); S. K Das, J. Alam, P. Mohanty 

and B. Sinha Phys. Rev. C 81, 044912 (2010); S. K. [29 

Das and J. Alam, arXiv: 1008.2643 [nucbthj; S. K. 

Das and J. Alam, arXiv: 1101.3385 [nucbth]. [30 
[13] Y. Akamatsu, T. Hatsuda and T. Hirano, Phys. Rev. 

C 79, 054907 (2009) [31 
[14] W. M. Alberico et al, Eur.Phys.J.C 71 1666 (2011) 
[15] S. Majumdar, T. Bhattacharyya, J. Alam and S. K. 



Das, Phys. Rev. C 84, 044901 (2011). 

J. Alam, S. Raha and B. Sinha, Phys. Rev. Lett. 

73, 1895 (1994). 

G. D. Moore and D. Teaney, Phys. Rev. C 71, 
064904 (2005). 

S. K Das, J. Alam and P. Mohanty, Phys. Rev. C 
82, 014908 (2010). 

M. Laine,JHEP 04 124, 201 1 . larXiv: 1 103.0372 1 [hep- 
ph]. 

M. He, R. J. Fries and R. Rapp, Phys. Let. B701, 
445 (2002) 

S. Ghosh, S. K. Das, S. Sarkar and J. Alam, Phys. 
Rev. D(R) 84,011503(R) (2011) 
L. Abreu, D. Cabrera, F. J. Llanes-Estrada and J. 
M. Torres-Rincon, Ann. Phys. 326, 2737 (2011). 
M. Djordjevic, M. Gyulassy, R. Vogt, S. Wicks, 

Nucl. Phys. A 774, 689 (2006). 

M. He, R. J. Fries and R. Rapp, larXiv:1106.6006l 
[nucl-th] 

P. Braun-Munzinger, K. Redlich and J. Stachel in 
Quark Gluon Plasma, R. C. Hwa (ed.) 491, e-Print: 
nucl-th/0304013 

B. Svetitsky, Phys. Rev. D 37, 2484( 1988). 

L.S. Ceng, N. Kaiser, J.Martin-Camalich and 

W.Weise, Phys. Rev. D 82, 054022 (2010) 

A. V. Manohar and M. B. Wise, Heavy quark 

physics, Camb. Monogr. Part. Phys. Nucl. Phys. 

Cosmol. 10, 1 (2000). 

Y.R. Liu, X. Liu, and S.L.Zhu Phys. Rev. D 79, 
094026 (2009) 

F.K. Guo, C. Hanhart, U. G. Meissner, Eur. Phys. 
J. A 40 171 (2009) 

H. Ohiki, H. Matsufuru, and T. Onogi, Phys. Rev. 
D 77, 094509 (2008) 



